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ABSTRACT 

All globular clusters (GCs) studied to date show evidence for internal (star-to-star) variation in their light 
element abundances (including Li, C, N, O, F, Na, Mg, Al, and probably He). These variations have been 
interpreted as evidence for multiple star formation episodes within GCs, with secondary episodes fueled, at 
least in part, by the ejecta of asymptotic giant branch (AGB) stars from a first generation of stars. A major 
puzzle emerging from this otherwise plausible scenario is that the fraction of stars associated with the second 
episode of star formation is observed to be much larger than expected for a standard IMF. The present work 
investigates this tension by modeling the observed anti-correlation between [Na/Fe] and [O/Fe] for 20 Galactic 
GCs. If the abundance pattern of the retained AGB ejecta does not depend on GC mass at fixed [Fe/H], 
then a strong correlation is found between the fraction of current GC stellar mass comprised of pure AGB 
ejecta, f p , and GC mass. This fraction varies from 0.20 at low masses (10 4 5 M Q ) to 0.45 at high masses 
(10 6 5 Mo). The fraction of mass associated with pure AGB ejecta is directly related to the total mass of the 
cluster at birth; the ratio between the initial and present mass in stars can therefore be derived. Assuming a star 
formation efficiency of 50%, the observed Na-O anti-correlations imply that GCs were at least 10-20 times 
more massive at birth, a conclusion that is in qualitative agreement with previous work. These factors are lower 
limits because any mass-loss mechanism that removes first and second generation stars equally will leave f p 
unchanged. The mass-dependence of f p probably arises because lower mass GCs are unable to retain all of the 
AGB ejecta from the first stellar generation. Recent observations of elemental abundances in intermediate-age 
LMC clusters are re-interpreted and shown to be consistent with this basic scenario. The small scatter in f p 
at fixed GC mass argues strongly that the process responsible for the large mass loss is internal to GCs. A 
satisfactory explanation of these trends is currently lacking. 
Subject headings: Galaxy: globular clusters — globular clusters: general 



1. INTRODUCTION 

Evidence has been accumulating for the past thirty years 
that globular clusters (GCs) harbor internal (star-to-star) vari- 
ation in their light e lement abundances (including Li, C, N, O, 
F. Na Mg. an d Al:|Cohenlll978UKrafdll979USrmth & Norris l 
1982b, 19 8"! lKraftlll994HGratton et al.ll2004t iPasquini et all 
2005; Smi th et alJ 120051) . And yet, save the most massive 
GCs, they can still be considered mon o-metallic in heavier 
elemen ts including Ca, Si, and Fe (e.g.. lCarretta et afl 2009a. 
l2010bl) . 

While early work focused on giant branch stars, abundance 
variations have now been observed in main sequence turn-off 
stars ( e.g.JGratton et aLll2001l:lBrilev et al.ll2002tlCohen et al.1 
12001 ICannon et al.[Tl998t iPancino et al l 1201 Obi) , and so the 
observed variations cannot be attributed to non-canonical 
mixing in evolved stars[]. Rather, the observed star-to-star 
variation in elemental abundances must be due to the fact that 
the stars formed from different material. 

Perhaps the most striking result emerging from these ob- 
servations is that the number of stars within a GC that show 
anomalous abundance ratios is comparable to those stars that 
have normal ratios. In this context 'normal' abundance ratios 
refers to abundances characteristic of field stars at the same 
[Fe/H] abundance, and 'anomalous' refers to abundance ratios 

1 Although non-canonical mixing cannot explain the totality of the ob- 
served star-to-star variation within GCs, there are well-documented corre- 
lations between some elemental abundances , such as c arbon and nitrogen, 
and location along the giant branch (see e.g., Smith 2002) that cannot be ex- 
plained by state-of-the-art stellar evolutionary models. This fact significantly 
complicates interpretation of certain elements. 



that differ markedly from the field. The comparable number 
of normal and anomalous stars is observed for all GCs studied 
to date , spanning a wide range in stellar mass and metallicity 
(e.g.. lMartell & Smithll2009tlCarretta et aTll2010d) . 

The anomalous stars display abundance ratios that pro- 
vide important clues to the source of the raw material from 
which they formed. These stars show enhanced Na and Al 
and depleted O and Mg abundances. They are also CN- 
enhanced and CH-depleted. These peculiar abundance pat- 
terns arise naturally when matter is brought to very high tem- 
peratures (far exceeding 10 7 K). At sufficiently high tem- 
peratures the CNO, Na-Ne, and Mg-Al nuclear reaction cy- 
cles are activated (the precise temperature required for acti- 
vation of these cycles depends in detail on the site, whether 
e.g., the stellar interior or the at the base of the convective 
envel o pe :ICharbonnel & Prantzosl2006tlKarakas & Lattanziol 
I2007t IVentura & D'Antonal l2008bl) . In intermediate-mass 
asymptotic giant branch (AGB) stars, the base of the convec- 
tive envelope reaches temperatures that are necessary to ig- 
nite these nuclear cycles. Because the envelope is convective, 
material through the whole envelope can therefore participate 
in nuclear burning. This process is known as envelope, or 
hot bottom burning, and it occurs for stars of initial masses 
in the range 4 M m < M < 8 M W (e.g.. lRenzini & Volil[l98Tt 
IVentura & D'An tona 2008b). Formation of the anomalous 
abundance ratios within the envelopes of AGB stars is appeal- 
ing because these stars do not produce heavier elements such 
as Ca, Si, and Fe, which then explains the lack of observed 
variation in these elements in most GCs studied to date. These 
stars also produce large amounts of He, which seems neces- 
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sary to explain the CMD morphology of many GCs. 

While AGB stars are plausible candidates for the source 
of the peculiar abundance patterns, other sites have been 
proposed. Alternatives include massive (> 20 Mq) rotat- 
ing stars dDecressin et al.l [2007b!) and massive binary stars 
(de M ink et al.ll2009l) . One of the key differences between the 
massive star and AGB scenarios is the timescales involved. In 
the massive star scenario the peculiar abundance patterns are 
created on a timescale comparable to the production of type 
II supernovae (SNe). This scenario therefore faces two addi- 
tional difficulties that the AGB scenario does not: 1) how to 
retain the processed material within the natal GC in the face 
of energy injection from SNe; and 2) how to create peculiar 
abundance patterns that do not show evidence for type II SNe 
products such as Si an d Fe. Other arguments aga inst this sce- 
nario are discussed in Conroy & Spergel (1201 ll) . While the 
AGB scenario is currently the favored source of the peculiar 
abundance patterns, it should be noted that this scenario con- 
tains significant uncertainties and theoretical difficulties. 

In light of these considerations, an intricate scenario has 
emerged to explain the observations. The key ingredient 
is that multiple generations of star formation have occurred 
within GCs. Later generations of stars formed from gas en- 
riched by the ejecta of AGB stars associated with earlier 
generations of star formation. Such a sce nario has been re- 
cently reviewed and discussed at length bv lConrov & Spergel 
d2011[). to which the read er is referred for d etails (see also, 
e.g.. [Cottrell & Da Cost J 1198 U lSmitj3ll987t i D'Ercole etail 
2008; iRenzinil 120081; iCarretta et all 1201 Oct iD'Ercole et ail 
2010). As discussed in that work, the number of star for- 



mation events in typical GCs was probably limited to two 
— a first generation of stars that formed from material with 
abundance patterns similar to field stars, and a second gener- 
ation formed from gas enriched by AGB ejecta from the first 
generation. The timescale between these two epochs of star 
formation is probably several 10 8 yr, which is the timescale 
for several relevant physical processes including AGB evo- 
lution, the time it takes for UV photon production to drop 
enough to allow the gas to cool and form stars, and the on- 
set of prompt type la supernovae. This small difference in 
age cannot be seen in the main sequence turn-off point in the 
ancient GCs, but can be observed in younger clusters. In- 
deed, the massive intermediate-age clusters in the Large Mag- 
ellanic Cloud (LMC) have a spread in their turn-off point that 
is consistent within an i nternal age spread of several 10 s yr 
dGoudfrooii et all 120091; iMilone et all 120091) . suggesting that 
this scenario occurs at the present epoch as well as the distant 
past. 

Another important ingredient in this scenario is that 
pure AGB ejecta must be mixed with gaseous mate- 
rial that has abundance ratios similar to the first genera- 
tion dPrantzos et al.ll2007t I Ventura & D'Antonall2008al, 120091; 
ID'Ercole et al.ll2010l) . This is required to reproduce the ob- 
served range in light element abundances, which extends 
smoothly from the very anomalous to the normal. A natu- 
ral way to acquire additional material with normal abundance 
ratios i s via accretion from the ambi e nt inte r stellar medium 
(ISM: iPflamm-Alte nbure & Kroupal 120091; ID'Ercole et aT] 
120 10t [Conrov & Spergel 201 1 ), although other scenario s are 
possible dD'Ercole et alj|2008t ; iGrat ton & C arrettall20lol) . As 
demonstrated in IConroy & Spergel 120 1 ll) . significant accre- 
tion from the ambient ISM is possible for the physical condi- 
tions characteristic of young GCs (i.e., cold dense interstellar 



media and low relative velocities). Ram pressure is not im- 
portant except for the lowest masses, where indeed anoma- 
lous abundance ratios are not observed, both in the Galac- 
tic open clusters dde Silva et alJl2009t iMartell & Smithll2009l: 
iPancino et al.ll2010ah and the in termediate-age clusters in the 
LMC (Conroy & Spe rgel 2011). Implicit in this requirement 
is that the accreted material remain incompletely mixed with 
the AGB ejecta — the stars must form from material with a 
range of abundances, not just some average of accreted mate- 
rial and AGB ejecta. 

The greatest challenge facing this otherwise plausible sce- 
nario is in explaining the roughly comparable number of first 
and second generation stars, which is observed for GCs span- 
ning a wide range in mass. Under the assumption that sec- 
ond generation stars form from AGB ejecta plus a modest 
amount of accreted material, the ratio of first to second gener- 
ation stars should be of order 10:1. That is, the standard sce- 
nario predicts a number of second generation stars lower by a 
factor of ten compared to observations. The standard predic- 
tion assumes 100% star formation efficiency and is based on 
stellar evolution theory and a canonical initial mass function 
(IMF), which implies that only ~ 10% of the mass of a stellar 
population ends up in AGB ejecta. A variety of solutions to 
this problem have been proposed, including a different IMF 
between first and second generation stars dSmith & Norrisl 
Il982at iPrantzos & Charbonnell 120061). and a substantially 
larger mass at birth of all GCs j5 'Antona & Caloi 2004; 
Bekki & Norrisl 120061. ID' Antona et alJl2007t iDecressin et all 
2007at ID'Ercole et all 120081; ISchaerer & Charbonnell 1201 lh . 



The second scenario requires that GCs were factors of 10- 
100 more massive at birth. If this is correct, it constitutes a 
dramatic revision of our understanding of GC formation and 
evolution. This tension provides the motivation for the present 
analysis. 

In the present work it is assumed that AGB stars with 
masses in the range 3-8 Mq contribute to the formation of 
the second generation of stars. However, the current genera- 
tion of AGB nucleosynthetic yields favors a narrower range of 
stellar masses, perhaps only 5-8 Mq, that are contributing to 
form the observed range of abundance patterns of second gen- 
eratio n stars (e.g.. lVentura & D'Antonall2009t ID'Ercole et all 
2010). Appealing to a smaller range of masses only exacer- 
bates the problem noted above because an even larger initial 
GC is required to produce the same amount of polluted ma- 
terial (a factor of two larger if one considers only 5-8 Mq 
AGB stars, rather than 3-8 M© AGB stars, everything else 
being equal). Moreover, if only massive AGB stars are al- 
lowed to contribute, then the timescale for the formation of 
second generation stars shrinks to < 10 8 yr. 

If a second generation of stars forms within the poten- 
tial well of a first generation GC, one might expect the 
second generation stars to be more spatially concentrated 
than the first generation, at least initially. A^-body simula- 
tions have shown that if they begin more centrally concen- 
trated, second generation stars will remain more concentrated 
than the first generat i on for a few central relaxation times 
dDecressin et al.l2008l: ID'Ercole et alj2008l) . after which time 
the relative number of first and second generation stars be- 
comes constant within the half-mass radius. There is evi- 
dence that the second generation is more centrally c oncen- 
trated than the first in NGC 1851 dZoccali et al.l l2009), wCen 
dPancino et al.l 12001 ISoTiima et all I2007l)~and NGC 3201 
(Carretta et al. 2010 a]). Th e situation for other GCs is less 
clear (ICarretta et al.ll2009d) . owing to the small numbers of 
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stars observed, incomplete spatial coverage, and the short re- 
laxation times for many GCs. 

Parallel to the advances in the chemical abundances of nor- 
mal GCs has been the revelation of distinct stellar popula- 
tions in the color-magnitude diagrams (CMDs) of the most 
massive GC s as revealed by the Hubble Space Telescope (see 
iPiottol 120091 for a review). It is now clear that the massive 
GCs uiCen and NGC 2808 harbor multiple distinct main se- 
quences as seen in their CMDs. Many more GCs harbor dis- 
tinct sub-giant branches and significant width in their red gi- 
ant branches. Although it is clear that the spread in the CMDs 
is related to the elemental abundance variations dYong et alJ 
20081: iMarino et alj|2008t ICarretta et al.1l2009ct iMilone et al.1 



20101) . a comprehensive understanding of the relation between 



these two phenomena is currently lacking. In particular, while 
much attention has been focused on the peculiar CMDs of the 
most massive GCs, it is far from clear that the features of these 
GCs are characteristic of the population as a whole. 

The present work aims to provide a quantitative understand- 
ing of the relative frequencies of first and second generation 
stars, the amount of material accreted from the ambient ISM, 
and the amount of AGB ejecta required to explain the ob- 
served elemental abundance variations, over a factor of 100 
in GC mass. This work therefore aims to bridge our under- 
standing of the properties of lower mass clusters with the most 
peculiar massive GCs. One of the major goals of the present 
analysis is to answer the following question: "how much more 
massive did the progenitors of present day GCs need to be 
in order to explain the observed abundance variations within 
GCs?" 

2. A SIMPLE MODEL FOR THE ABUNDANCE VARIATIONS 

One of the most intriguing observational results is the 
strong anti-correlation between [Na/Fe] and [O/Fe] within 
a given cluster. It has been shown in previous work 
that this anti-correlation can be reproduced with a sim- 
ple dilution model wherein second generation stars form 
from nor mal and processed material mixed in varying 
amounts dPrantzos et all l2007t IVentura & D'Antonal 2008a, 
l2009t iD'Ercole et al.ll2010i) . Processed material means in this 
context matter that has been exposed to temperatures exceed- 
ing «7x IQ 7 K, where the CNO and Na-Ne nuclear reactions 
are active in the convective envelopes of massive AGB stars. 
Stars formed from pure processed material have enhanced 
[Na/Fe] and depleted [O/Fe] abundances. Normal material 
means here matter that has the same abundance patterns as 
the initial stellar population, and also probably has the same 
abundance patterns as field stars in the Galaxy, at the same 
[Fe/H]. Normal material at the low metallicities considered 
has enhanced [O/Fe] and either solar or sub-solar [Na/Fe] 
abundances. In the present discussion, AGB ejecta' will be 
synonymous with 'pure processed material', as it is assumed 
herein that the source of the processed material is AGB ejecta. 

If the normal and processed abundances are denoted with 
subscripts 'o' and 'p' and if the fraction of pure processed 
material in the j— th star is f p then the abundance of element 
i in a star formed from a mixture of normal and processed 
material is simply: 



[i/¥c]j = log ( (1 -ft)lQV/**+fJlQV/*>h 



(1) 



assuming Fe is constant. 
It is important to notice that in the context of this model 



there will be stars with f p = 0.0; these are truly first gener- 
ation stars that formed purely from normal material. Subse- 
quent generations of stars then form from a mixture of AGB 
ejecta and normal material. The source of the normal material 
required to produce the observed [Na/Fe] and [O/Fe] values 
in the second generation is currently unknown. A plausible 
source is accretion from the ambient ISM over the several 10 8 
yr when intermediate-mass AGB stars are evolving, and af- 
ter the type II SNe h ave exploded and cleared the GC of any 
remaining initial gas (IConrov & Sp ergel 201j]). 

Once the normal and processed abundances of [Na/Fe] and 
[O/Fe] are specified, one can readily estimate fj, for each star 
in a GC by comparing the star's observed Na and O abun- 
dances to the model predictions. An estimate of the global 
value of f p can then be obtained by averaging the individ- 
ual // values over all the stars in the GC. This global f p rep- 
resents the fraction of the present stellar mass comprised of 
pure AGB ejecta. The fraction of the present GC stellar mass 
formed from normal material is then given by 1 —f p . 

This procedure for determining f p is straightforward. An 
advantage of this approach is that no arbitrary, sharp distinc- 
tion is made between first and later generations of stellar gen- 
erations based on their location in the Na-0 plane. This is 
advantageous because the observed distribution of stars along 
the Na-O anti-correlation is continuous in most GCs, and so 
there is no natural way to separate first and later generations. 

A complication in interpreting the results from such a dilu- 
tion model is that the fraction of normal material is not eas- 
ily separable into the normal material comprising the initial 
stellar population and the normal material locked in second 
generation stars. This ambiguity arises because stars with 
< f p < 0.1, which are second generation stars, are con- 
fused with stars with f p = 0.0, which are truly first genera- 
tion stars. The fraction of normal material that is associated 
with second generation stars will be denoted by / acc . The 
fraction of total mass in the first stellar generation is then 
fx = (1 — fp){\ — /acc)> an d the fraction of total mass in accreted 
material is f a = (1 —f p )f aC c. The fraction of mass in the second 
generation is then = f a +f p . 

Ultimately we are interested in knowing if the amount of 
mass in pure AGB ejecta incorporated in second generation 
stars and given by Mf p where M is the total GC stellar mass, 
can be produced from the first stellar generation. The initial 
mass needed to produce the observed amount of AGB ejecta 
can be estimated as follows. First, define /agb as the frac- 
tion of initial mass that ends up in AGB ejecta and is avail- 
able for second generation star formation. It can be computed 
from stellar evolution in conjunction with an IMF. Next, de- 
fine 6sf as the star formation efficiency, or the fraction of gas 
mass within the young GC that ends up in second generation 
stars. The initial mass then needed to produce the observed 
amount of AGB ejecta is M/p//agb/£sf- The ratio between 
this quantity and the actual amount of mass in long-livecQ first 
generation stars, M(l —f p )(l -/acc), is: 



/mi = max[ 1.0, 



fp 



£Sf/aGb(1 -/p)(l "/acc) 



(2) 



2 For a canonical Krourja (2001) IMF, approximately 50% of the initial 
stellar mass of a population is lost due to stellar evolution and death. In the 
present analysis attention is focused on long-lived stars, so the massive stars 
that carry away substantial mass when they die are not included in bookkeep- 
ing except where explicitly mentioned. 
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Table 1 

Summary of Symbols and Quantities 



Derived From Description 



The quantity, /mi is the mass enhancement factor for the 
first generation, considering only long-lived stars. If stars of 
all masses are included then /mi should be increased by a fac- 
tor of w 2 (see footnote 2). The definition of /mi in Equa- 
tion |2]preserves the requirement that /m\ cannot be less than 
unity. It is instructive to consider a quantitative example. If 
/agb ~0-l ( as suggested by stellar evolution and a canonical 
iKroupal J200l IMF), f p = 0.5, / acc = 0.5, and e SF = 0.5, then 
/mi =40.0. This means that the first generation contained 40 x 
more long-lived stars at birth compared to its present mass. 

In what follows a star formation efficiency of esF = 0.5 will 
be adopted. The simplest argument for a high star formation 
efficiency in GCs is based on the Virial Theorem and the re- 
quirement that GCs remain gravitationally bound. This ar- 
gument suggests that the total mass in stars formed needs to 
be at least 50% of the initial gas mass (the resulting bound 
mass depends on h ow quickly gas is removed from the sys- 
tem, for details see lLada et al] 1 19841) . The gas available for 
second generation star formation finds itself in a unique con- 
figuration in that it is at the bottom of a deep potential well 
(set by the first generation stars). In principle the star forma- 
tion efficiency could therefore be lower and yet still result in 
a bound cluster. The uncertainty carried by the parameter 6sf 
is an inevitability of the present analysis. 

While recognizing that in reality star formation is not 100% 
efficient, previous work has nonetheless made the assumption 
that 100% of AGB ejecta from the first generation is incorpo- 
rated into second generation stars. This assumption was made 
for the simple reason that it minimizes the required mass en- 
hancement factor. In the present work a more realistic value 
for £sf has been adopted. 

Another quantity of interest is the total mass enhancement 
factor in long-lived stars, that is, the ratio of total mass at birth 
to total present mass. This quantity is 

f t =mn(l.0,—&—\ (3) 

V ESF/ AGB / 

Notice that this is smaller than f m by a factor of (1 — f p )(l - 
/acc) for / > 1. In the example above, the total mass enhance- 
ment factor would be 10, a factor of four smaller than the 
mass enhancement factor for the first generation. The former 
is smaller than the latter because of the addition of second 
generation stars, which at least partially compensates for the 
substantial loss of first generation stars. For reference, a GC 
evolving in isolation that experiences no secondary genera- 
tions of star formation and also experiences no loss of stars 
from e.g., stellar evaporation or ejection, will have /mi = 1 
and / = 1. 



Most authors assume that second generatio n stars are born 
with masses in the range 0.1-0.8 M Q (e.g. jD'Ercole et all 
120081: iLind et all 1201 ll) . This assumption is made in order 
: to minimize the required mass enhancement factors. If AGB 
.ejecta is only placed in long-lived second generation stars, 
then none of the ejecta is 'wasted' on higher mass stars that 
evolve and die on short timescales. This assumption reduces 
the required mass enhancement by a factor of ps 2 compared 
to the assumption, made in the present work, that seco nd gen- 
eration stars fully populate a standard Kroupa (2001) IMF. A 
fully populated IMF is assumed here because there is no phys- 
ical reason to believe that the IMF should be biased toward 
low masses during second generation star formation. 

An argument often made for an IMF of second generation 
stars that is skewed toward low masses is that massive second 
generation stars would explode and thereby truncate further 
star formation. This is certainly an issue that must be ad- 
dressed in any comprehensive theory, but i t is not a strong ar- 
gumen t for a skewed IMF. As discussed in lConrov & Spergell 
d201 ll) . the heating rate of UV photons should be substantial 
during the first ~ 10 8 yr of the young GC. This heating may 
be sufficient to delay star formation until the UV photon pro- 
duction rate drops precipitously at ~ 10 8 yr. After this time, 
the gas may cool catastrophically in a manner analogous to 
the formation of the first generation of GC stars. This simple 
scenario suggests that massive second generation stars could 
form and yet not have a debilitating effect on the conversion 
of a significant fraction of the accumulated gas into stars. 

The symbols and quantities defined in this section are sum- 
marized in Table[T]for convenience. 

3. CONSTRAINTS ON THE BIRTH MASSES OF GALACTIC GCS 

Recently, Carretta and collaborators have published abun- 
dance measurements for > 1800 red giant stars in 20 GCs 
spanning a wide range in mass and metallicity dCarrettal 
l2006t ICarretta et af]|2007al I2009d I2010dl) . The abundances 
were derived from high resolution spectra obtained with the 
FLAMES spectrograph at the VLT. Basic data for this sample 
are collected in Table [2] including common names, average 
[Fe/H] values, and number of stars with abundance determi- 
nations or upper limits for both [O/Fe] and [Na/Fe]. Stellar 
masse s are derived from luminosities adopted from the lHarrisl 
( 1996) catalog assuming M/Ly = 2 M Q /£.©. The table also in- 
cludes a variety of derived products discussed in detail below, 
including the inferred [Na/Fe] and [O/Fe] abundances of both 
processed and normal material, the fraction of processed ma- 
terial (AGB ejecta), the average distance of the stars observed 
to the cluster center in units of the half-mass radius (the latter 
adopted from the Harris catalog), and the mass enhancement 
factors. 

Stars were used in the following analysis only if they have 
reported abundance measurements or upper limits for both 
[Na/Fe] and [O/Fe]. Some stars were observed in a configura- 
tion that did not include the spectral region covering the [OI] 
line; such stars are not included in the analysis below. The 
decision of whether or not to include the [OI] line is not cor- 
related with the actual [O/Fe] abundance, so removing such 
objects should not affect the results. 

The model [Na/Fe] and [O/Fe] abundances of the normal 
and processed material for each cluster were estimated by 
eye from the distribution of stars in the [Na/Fe] -[O/Fe] plane. 
The estimated abundances of the normal material lie within 
the range of the abundances of field sta rs at the same [Fe/H] 
dVenn et al.ll2004t ICarretta et al.ll2010cT) . 



fp 


Na-O dilution model 


fraction of GC mass comprised of AGB ejecta 


face 


free parameter 


fraction of normal material accreted 






from ambient ISM 


/agb 


stellar evolution, IMF 


fraction of mass in AGB ejecta 






from a coeval stellar population 


c Sf 


free parameter 


efficiency of star formation 


fa 


(!-/»)/«« 


fraction of GC mass in accreted material 


f\ 


(1-/ P )(l-Acc) 


fraction of GC mass in first generation stars 


h 


/«.+/? = 1-/1 


fraction of GC mass in second generation stars 


.hi] 


EquationE] 


mass enhancement factor for first generation 


ft 


Equation]!] 


total mass enhancement factor 
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Table 2 

Summary of GC Data 



ID 




[Fe/H] 


N 


log(M) 


[Na/Fe] 


[0/Fe] 


[Na/Fe] p 


[0/Fe] p 


fp 


(R)/R h 


Jm ] 


/, 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


NGC 7099 


M30 


-2.33 


19 


5.19 


-0.30 


0.50 


0.80 


-0.80 


0.32 


2.33 


18.7 


6.37 


NGC 7078 


M15 


-2.33 


20 


5.89 


-0.30 


0.50 


0.80 


-0.80 


0.36 


1.47 


22.4 


7.19 


NGC 4590 


M68 


-2.23 


36 


5.16 


-0.35 


0.60 


0.80 


-0.80 


0.28 


1.83 


15.2 


5.52 


NGC 6397 




-1.98 


13 


4.87 


-0.30 


0.40 


0.80 


-0.80 


0.20 


1.04 


10.1 


4.04 


NGC 6809 


M55 


-1.98 


75 


5.24 


-0.25 


0.45 


0.80 


-0.80 


0.33 


1.14 


19.9 


6.65 


NGC 6715 


M54 


-1.57 


76 


6.23 


-0.15 


0.40 


0.80 


-0.80 


0.42 


3.65 


29.5 


8.49 


NGC 1904 


M79 


-1.55 


39 


5.37 


-0.25 


0.30 


0.80 


-0.80 


0.31 


1.92 


17.9 


6.18 


NGC 6752 




-1.56 


■SS 


5.31 


-0.10 


0.50 


0.80 


-0.80 


0.36 


1.69 


22.1 


7.12 


NGC 6254 


M10 


-1.56 


77 


5.21 


-0.40 


0.40 


0.70 


-0.80 


0.27 


1.54 


14.9 


5.44 


NGC 3201 




-1.50 


94 


5.21 


-0.30 


0.30 


0.60 


-0.80 


0.34 


1.12 


20.2 


6.72 


NGC 5904 


M5 


-1.34 


106 


5.75 


-0.20 


0.45 


0.70 


-0.80 


0.38 


1.73 


24.5 


7.59 


NGC 6218 


M12 


-1.31 


66 


5.15 


-0.20 


0.50 


0.80 


-0.80 


0.34 


1.62 


20.9 


6.87 


NGC 288 




-1.23 


64 


4.92 


-0.10 


0.25 


0.80 


-0.80 


0.29 


1.35 


16.5 


5.83 


NGC 6121 


M4 


-1.20 


80 


5.10 


-0.10 


0.40 


0.80 


-0.80 


0.33 


0.88 


19.3 


6.52 


NGC 6171 


Ml 07 


-1.06 


27 


5.07 


0.00 


0.40 


0.80 


-0.80 


0.31 


2.01 


17.7 


6.14 


NGC 2808 




-1.10 


90 


5.98 


-0.10 


0.40 


0.60 


-0.80 


0.42 


3.14 


29.0 


8.40 


NGC 6838 


M71 


-0.80 


31 


4.46 


0.00 


0.50 


0.80 


-0.80 


0.25 


0.93 


13.3 


4.98 


NGC 104 


47Tuc 


-0.74 


109 


5.99 


0.10 


0.35 


0.75 


-0.80 


0.40 


1.66 


27.2 


8.10 


NGC 6388 




-0.40 


29 


5.99 


0.00 


0.20 


0.70 


-0.80 


0.39 


5.56 


26.1 


7.90 


NGC 6441 




-0.34 


24 


6.08 


-0.10 


0.20 


0.70 


-0.80 


0.36 


9.16 


23.0 


7.30 



Note. — (1): average GC [Fe/H] abundance; (2): number of stars used in the analysis; (3): logarithm of GC stellar mass in 
units of Mq; (4): [Na/Fe] abundance of normal material; (5): [O/Fe] abundance of normal material; (6): [Na/Fe] abundance of 
pure processed material; (7): [O/Fe] abundance of pure processed material; (8): AGB ejecta mass fraction; (9): average cluster- 
centric radius of stars in units of the cluster half-mass radius; (10): first generation mass enhancement factor; (11): total mass 
enhancement factor. 



The distribution of [Na/Fe] and [O/Fe] abundances is shown 
for all 20 GCs in Figure [Tj along with the corresponding di- 
lution models. The typical error on the [Na/F e] and [O/Fe] 
abunda nces is 0.08 and 0.14 dex, respectively (Carret ta et al.1 
2009c). These errors are dominated by uncertainties in the 
measured equivalent widths, but also include error due to 
the adopted atmospheric parameters. The GCs are sorted by 
increasing stellar mass. It is apparent that the lower mass 
GCs have on average shorter Na-O anti-correlations than the 
high er mass clusters, a trend tha t has been noticed before 
(e.g. JCarretta et al1l2007bl [201 Od) . This trend, readily appar- 
ent in the data, will be the ultimate source of the GC mass- 
dependent trends discussed below. 

3.1. Derivation of f p for Galactic GCs 
By comparing the abundance patterns of stars within each 
GC to the corresponding dilution model, f p for each star and 
hence the average f p for the cluster can be estimatecQ. The 
result is shown as a function of GC stellar mass in Figure 
[2] The parameter f p is only weakly sensitive to the adopted 
abundances of the pure processed and normal material in the 
dilution model, as discussed below. Errors on f p have been 
estimated in the following way. New dilution models were 
run for each GC with variation in [0/Fe] p by ±0.4, [Na/Fe] ;; 
by ±0.1, [0/Fe] by ±0.05, and [Na/Fe]„ by ±0.1, and the 
resulting minimum and maximum f p values from this proce- 
dure were adopted as the standard error. The range adopted 
for each parameter was determined by the typical variation 
allowed by the data. The [Na/Fe] ; , abundance is relatively 
well-constrained by the data, especially for GCs with stars 
both above and below the knee in the model. The choice of 
[0/Fe] p has only a modest effect on f p because of the nature 

3 For stars with upper limits on [O/Fe] abundances, only the [Na/Fe] values 
were used to determine fl . Other choices for how to handle upper limits do 
not change the overall results. 



of the dilution model - stars on the upper branch of the Na-O 
correlation can have very different [O/Fe] values and yet very 
similar //. 

There is a strong correlation between GC mass and the 
fraction of mass comprised of pure AGB ejecta. The cor- 
relation is particularly strong when restricted to GCs with 
[Fe/H] > -1 .5. This correlation is a quantitative manifestation 
of the trend noted above that low mass GCs have short Na- 
O anti-correlations compared to high mass GCs. The best-fit 
relation between / ; , and GC mass is: 

f p = (0.10±0.021)log(M/M Q )-(0.21 ±0.12). (4) 

There are only two GCs at the low mass end, and one might 
wonder to what extent they are influencing the observed cor- 
relation. If those two clusters are removed, the significance 
of the correlation is reduced from 4.8cr to 3.1er; the trend is 
therefore robust to the removal of the two lowest mass clusters 
in the sample. 

As discussed in the Introduction, there is some evidence 
that the extent of the Na-O anti-correlation is a function of 
radius within GCs. One might then wonder if the trend seen 
in Figure [2] is due to a selection effect. In order to test for 
such an effect, for each GC the average radius of stars for 
which abundance data are available was computed in units of 
the GC half-mass radius. The results are tabulated in Table [2] 
and shown as a function of GC mass in Figure [3] The lowest- 
mass GC and the four GCs with (R)/Rh > 3 define a weak 
trend with mass. If these five GCs are excluded, there is no 
remaining trend with mass in Figure[3] and yet the strong cor- 
relation between f p and mass remains. Future observations 
over a greater range in radius would be valuable, but even 
with current data it is clear that the trend of f p with mass is 
not due to sampling systematically different regions of GCs 
as a function of their mass. 

As can be seen in Table [2] the adopted pure AGB ejecta 
abundances [0/Fe] p and [Na/Fe] /; vary only weakly, or not at 
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[O/Fe] 

Figure 1. Na-O anti-correlation for 20 GCs. Data are from lCarrettj <2006]) , lGratton et alj j200l) . and lCarretta eTaTI )2007dlcll2009clbll2010j) . GCs are sorted 
by mass, and the logarithm of the GC stellar mass in solar units is shown in each legend, as is the average [Fe/H] of each cluster. Arrows indicate upper limits on 
[O/Fe] abundances. A typical error on the abundances is shown in the upper left panel. Lines show the dilution models, and open squares mark the location at 
which the contribution from AGB ejecta and normal material is equal. 
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Figure 2. Relation between the fraction of GC mass comprised of AGB 
ejecta, and GC stellar mass. GCs are color-coded according to their metal- 
licity. The solid line is the best-fit linear relation. The metal-rich bulge GC 
NGC 6441 is labeled. 

all, with GC mass and [Fe/H]. The abundances were chosen 
to be [0/Fe] p = -0.8 and [Na/Fe] p = 0.8 initially, and only 
changed if demanded by the data. The adopted abundances 
for each GC are not unique, and for the GCs with the shortest 
Na-O anti-correlations, there is a substantial amount of room 
for variation. 

The pure processed abundances coul d have been adopted 
from the latest AGB yields of IVentura & D'Antona 
(2009) , but, as noted b y those authors and others (e.g., 
Ventura & Marigo 2010), the yields are quite uncertain, 
especially for Na and O, and do not provide good fits to 
the observed a bundances of anomalous stars. For example, 
iD'Ercole et al.l (120121) were able to model the very low [O/Fe] 
abundances in NGC 2808 only after invoking an additional 
deep-mixing mechanism in red giants to lower the [O/Fe] 
abundances below the standard AGB yields. It is clear that 
AGB models still cannot, from first principles, produce 
elemental yields that match even the range of observed 
abundances in the second generation of stars in GCs. It is 
for this reason that the model AGB yields were not adopted 
when choosing [0/Fe] ;; and [Na/Fe] p . 

The pure AGB ejecta abundance ratios could also have been 
chosen for each GC to coincide with the most anomalous ob- 
served stars in each cluster. This would result in an increase 
in f p for each cluster. Consider for example NGC 6838, 
which as can be seen in Figure [T] contains a short Na-O anti- 
correlation. If the pure AGB ejecta abundances were chosen 
to be [O/Fe]^ = 0.0 and [Na/Fe]= 0.6, then the resulting pure 
processed fraction would be f p = 0.4 rather than / = 0.25. Do- 
ing this for each GC would yield a pure processed fraction, f p 
that was essentially independent of GC mass. 

However, there is no known reason why the processed 
yields should be a function of GC mass at fixed [Fe/H]. Com- 
pare for example NGC 2808 with NGC 6171. These GCs 
differ by a factor of ten in mass but have nearly the same 
[Fe/H] . The adopted pure processed abundances in NGC 6171 
can be justified based on the fact that the same abundances fit 



Figure 3. Average cluster-centric radius of stars for which abundance mea- 
surements are available, in units of the GC half-mass radius. Results are 
shown for all 20 GC as a function of GC mass. There is no strong trend with 
GC mass, which implies that the results in Figure [2] are not biased by ob- 
serving clusters of different masses at different average radii. The metal-rich 
bulge GC NGC 6441 is labeled as it is a significant outlier. 

the whole extent of the Na-O anti-correlation in NGC 2808. 
Comparing other GC pairs with similar metallicities, such as 
NGC 6838 and NGC 104, NGC 6397 and NGC 6809, or NGC 
6218 and NGC 5904, provides additional justification for the 
adopted pure processed yields. 

In summary, the derived pure AGB ejecta fractions, f p , and 
especially the trends with GC mass, are robust to the details 
of the model, but they do depend on the assumption that the 
AGB yields are independent of GC mass at fixed [Fe/H]. 

3.2. Characterization of the Multiple Stellar Population 
Phenomenon in Galactic GCs 

The relation between f p and GC mass derived in the previ- 
ous section can now be used to derive all of the other quan- 
tities described in Table Q] once the parameters esF, /acc and 
/agb are specified. S tandard stellar evolution in conjunction 
with a Kroupal (120011) IMF suggests /agb « 0.1 far AGB stars 
in the 3-8 M range, assuming that all of the AGB ejecta 
is later available for second generation star formation. Varia- 
tion of this parameter and / acc will be discussed below. A star 
formation efficiency of 6sf = 0.5 is also adopted. 

Relations between various derived properties and present 
GC mass are shown in Figure |4] Properties include the frac- 
tion of present GC mass in pure processed material, f p , in ac- 
creted material, f a , in first generation stars, f\, in the second 
stellar generation, fa, and mass enhancement factor of stars in 
the first generation, /mi , and total mass enhancement factor, 
ft- 

The left panel shows results assuming constant values for 
/agb and / acc . The right panel shows how the results change 
when these two parameters are allowed to be mass-dependent. 
The mass-dependence of these parameters was chosen so as to 
produce approximate mass-independence of the derived func- 
tions /mi, ft, fi, and f\. The variation in results between the 
left and right panels serves to illustrate the approximate range 
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Figure 4. Characterization of the multiple stellar population phenomenon in GCs. Properties shown include the fraction of present GC mass in pure AGB ejecta, 
fp, in accreted material, /„, in first generation stars, /[, in the second stellar generation, fi = fp + fa, fraction of second generation mass comprised of pristine 
accreted material, / acc , fraction of initial mass that ends up in AGB ejecta and is available for second generation star formation, /agb . mass enhancement factor 
in the first generation, f Mi , and total mass enhancement factor, ft . A star formation efficiency of egF = 0.5 is assumed. The thick black line is derived from the 
observations (see Figure|5J, the blue lines are assumptions, and the red and green lines follow directly from egp, /agb> /acc, and f p . In the upper portion of the 
figures, the y-axis on the left side describes the enhancement factors of long-lived stars, while the right side describes the enhancement including all stars. The 
approximate location of fyi for tdCen is indicated. Left panel: The parameters /agb and / aC c are assumed to be constant. Right panel: The mass-dependence of 
the parameters /agb ar| d /acc are chosen such that the majority of the other derived functions are independent of mass. 



in these parameters allowed by the data given the model un- 
certainties. 

It is entirely plausible that both /agb and / acc are in reality 
mass-dependent. For example, if Bondi accretion is the phys- 
ical process responsible for bringing in the accreted material, 
then / acc should increase with mass, while if simple geomet- 
ric cross sectional sweeping were domina nt, then it may de- 
creas e with increasing GC mass (see e.g.. lConrov & S pergel 
2011). The parameter /agb represents the fraction of stellar 
mass that is returned to the ISM as AGB ejecta and available 
for second generation star formation. If an increasing frac- 
tion of AGB ejecta is lost from the young cluster at lower 
GC masses, then /agb could decrease with decreasing mass. 
Such a scenario seems plausible since the wind velocity from 
AGB stars is in the rang e 10-20kms _1 dLoup et alj|1993b 
Vassiliadis & Wood 1993), which is of order the escape ve- 
locity for low mass GCs. 

The approximate value of f^i is shown for usCen, and is es- 
timated as follows. iRenzinil J2008) estimates that 0.7% of the 
initial mass of a stellar population would be returned to the 
ISM as fresh He (i.e., He produced within stars, as opposed to 
primordial He) from AGB stars with 3 < M < 8 M . In wCen, 
approximately 57% of the stars belong to the metal-poor com- 
ponent, 33% to the intermediate-metallicity component, and 



10% to the metal-rich component dPiotto et al.l 120051) . The 
CMD morphology of these components suggests He mass 
fractions of Y = 0.25 , Y ss 0.38, and Y ~ 0.4; the latter value 
is highly uncertain (ISollima et alJ l2005). Taken together, this 
implies a total mass in fresh He of sa 1.7 x 10 5 M Q in long- 
lived stars. If the metal-poor population represents the first 
generation, then it could have produced only 1.2 x 10 4 M Q 
of fresh He, assuming a total mass in long-lived stars of 
3x10 Mq for wCen. Assuming 6sf = 0.5, the mass enhance- 
ment factor for the first generation is therefore fu\ = 28. This 
independent estimate of /mi agrees fairly well with the value 
expected for its mass based on extrapolation of the model pre- 
sented herein. 

Several generic features stand out in Figure |4] First, the 
fraction of second generation stars is always > 50%, and i s 
only a weak function of mass (see also lCarretta et al.l l2010c). 
While fi is relatively constant, the composition of the sec- 
ond generation stars (given by f p and f a ) varies more strongly 
with mass, with the most massive GCs harboring second 
generation stars comprised largely of pure AGB ejecta, and 
lower mass GCs containing second generation stars com- 
prised mostly of normal material accreted from the ambient 
ISM. This is a natural consequence of the fact that lower mass 
GCs have shorter Na-0 anti-correlations compared to higher 
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mass GCs. 

Second, the mass enhancement factor for first generation 
stars is very large, approximately 30 for the most massive 
clusters, and at least 10 for the less massive ones. The to- 
tal mass enhancement factor, f t , is also quite large, though by 
factors of 2 - 3 smaller than /m\ ■ Previous work drew attention 
to the fact that fu\ has to be large in order to explain the ob- 
served abundance patterns, and therefore concluded that GCs 
must have been substantially larger at birth compared to their 
present masses. But this is only true with regards to the first 
generation population. Since the first generation is subdom- 
inant by number at the present epoch, the ratio of total mass 
at birth to the total present mass is factors of several smaller 
than /mi . 

In the present analysis attention has been focused on the 
mass enhancement factors for long-lived stars. However, 
short-lived stars (defined as stars with ages less than the 
present age of the Universe) comprise roughly one half of the 
t otal initial mass in a coeval stellar population formed from 
a lKroupal {2001) IMF. The mass enhancement factors for all 
stars are therefore a factor of approximately two larger than 
that quoted for long-lived stars only. For the most massive 
clusters, the total mass enhancement factors for all stars are 
therefore « 20. 

It is important to recognize that these mass enhancement 
factors are all strictly lower limits because any mass-loss 
mechanism that affects both first and second generation pop- 
ulations equally will not affect the distribution of stars in the 
Na-O plane, which is the fundamental observable underpin- 
ning the present discussion. For example, lower mass clus- 
ters (several 10 5 M Q ) can lose a significant fr action of their 
mass via stellar evaporation over a Hubble time dFall & Zhangl 
l200l . 

The basic conclusion from this section is that the ancient 
Galactic GCs had to be at least 10-20 times more massive at 
birth in order to produce enough AGB ejecta to account for 
the observed distribution of [Na/Fe] and [O/Fe] abundances. 

Finally, the ratio f p /f2 quantifies the fraction of mass in 
second generation stars comprised of pure AGB ejecta. This 
ratio is 30%-60% over the full range in GC mass. Broadly 
speaking, second generation stars are comprised of half pure 
AGB ejecta and half pristine material, in agreement with pre- 
vious modeling of the elem ental abundance variations within 
GCs dD'Ercole et al.ll2010h . 

4. ELEMENTAL ABUNDANCE VARIATIONS IN LMC CLUSTERS 

Knowledge of the abundance variations in star clusters out- 
side the Galaxy is limited because of the substantial distances 
to even our nearest neig hbors harboring clusters , the Mag- 
ellanic Clouds. Recently, Mucciar elli et al.l (120091) measured 
Na and O abundances for 18 giants in three old metal-poor 
GCs in the LMC. They found clear evidence for an Na-O anti- 
correlation in these clusters that looks broadly similar to the 
Na-O anti-co rrelations observe d in Ga lactic GCs. 

However, iMucciarelli et al.l (120081) measured elemental 
abundances in 27 red giant stars located within four 
intermediate-age (1-2 Gyr), moderate-metallicity LMC clus- 
ters, and concluded that the abundance patterns "show neg- 
ligible star-to-star scatter within each cluster". Recall that 
these and other intermediate-age LMC clusters show evidence 
in their main sequence turn-off points for an internal age 
sprea d of several 10 8 yr (Goudfrooij et al. 2009; Milon e et"afl 
2009). These results have been interpreted bv lBekkil (120111) 
as evidence for a qualitatively different scenario at work in 
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Figure 5. Na-O anti-correlatio ns for intermediate-age L MC clusters sorted 
by cluster mass. Data are from Mucciarelli et al. (2008). Stellar masses are 
estimated as described in the text and are shown in the legends in solar units. 
Dilution models are shown as lines. For each cluster, the data are subdivided 
by the mean sodium abundance. The mean abundance pattern for each sub- 
population is shown as stars. 

these clusters (wherein these clusters capture giant molecu- 
lar clouds that provide the fuel for secondary star formation, 
without AGB ejecta). It would be much more appealing, on 
the principle of simplicity, if the multiple stellar population 
phenomenon observed in these intermediate-age LMC clus- 
ters were similar to what is conjectured to have occurred in 
the ancient Galactic and LMC GCs. A thorough understand- 
ing of these clusters is also essential if one hopes to observe 
present-day young clusters in the process of forming second 
generation stars. 

Motivate d by these consid e ration s, the abundance data pre- 
sented in Mucciarel li et alj ((2008) are reinterpreted in the 
context of the dilution models presented herein. Figure [5] 
shows the Na-O anti-correlations for the four intermediate- 
age LMC clusters. Each cluster has an average [Fe/H] abun- 
dance of -0.3 to -0.5, which, for solar abundance ratios, corre- 
sponds to Z = 6 x 10~ 3 - 1 x 10~ 2 . The stellar masses of these 
clusters are est i mated by adopting V-band magnitudes from 
van den Berghl(fT98l . an average E(B — V) = 0.08, and a dis- 
tance modulus of 18.5. A V-band mass-to-light ratio of 0.4 
was adopted, which is appropriate for a stellar age of 1.6 Gyr 
andZ = 6xl0 -3 . 

The stars within each cluster have been split into two groups 
according to the mean sodium abundance of the stars. The 
mean sodium and oxygen abundances were then computed 
within each subpopulation. In all clusters there is statistically- 
significant evidence for different mean sodium abundances in 
the two populations. The two populations are most visually 
striking in NGC 2173 and NGC 1978. Of course, computing 
the mean abundance of subpopulations split according to that 
abundance will tend to exaggerate the differences between the 
two populations. Another way of assessing the statistical sig- 
nificance of the internal sodium variation is to consider the 
stars at the extremes of the distribution. In NGC 1783 the 
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stars with the highest and lowest [Na/Fe] abundances are con- 
sistent within la, while for NGC 1651, NGC 2173, and NGC 
1978 the most extreme stars differ at the > 3a level. There 
is therefore strong evidence for internal variation in [Na/Fe] 
in the clusters NGC 1651, NGC 2173, and NGC 1978, and 
ambiguous or weak evidence for NGC 1783. 

The [O/Fe] abundances show co nsiderably less va r iation , 
and it is perhaps for this reason that M ucciarelli et al.l (2008) 
concluded that there was negligible star-to-star variation 
within each cluster. However, two important facts conspire 
to reduce the expected variation in [O/Fe] in clusters that con- 
tain second generation star formation from AGB ejecta. First, 
the first generation stars in the intermediate-age clusters are 
not a-enhanced, and so the maximum value of [O/Fe] is re- 
duced to roughly 0.0. Second, at the moderate metallicities 
characteristic of these clusters, the depletion of O abundances 
due to hot bottom burning in AGB stars is much less than at 
lower metallicities characteristic of the ancient clusters. The 
mass- and metallicity-dependent A GB yields from the mod- 
els of[Ventura & D' Antonal d20 09) make this point clearly. In 
their Z = 4 x 10" 3 models, the [O/Fe] yields for AGB stars 
with masses in the range 3-6.5 M Q never drop below 0.0. 
So while the full range in [O/Fe] for the low metallicity GCs 
routinely exceeds 1 dex, at higher metallicities the range is 
expected to collapse to nearly zero, as observed. 

In Figure|5] schematic dilution models are included to guide 
the eye. Only NGC 1978 shows signs of a true Na-O anti- 
correlation. It is intriguing that this is also the most massive 
intermediate-age cluster in the sample. However, th is is the 
only massive cluster studied by M ilone et al.l (120091) where a 
significant spread in its main sequence turn-off point was not 
detected, suggesting a coeval population. As noted in that 
work, the CMD data for this cluster were obtained in the clus- 
ter outskirts, unlike the other clusters studied. The lack of an 
observed age spread could therefore be explained if the sec- 
ond generation is more centrally concentrated than the first 
generation, as models predict. 

While it is clear from the data that elemental abundance 
variations exist in the intermediate-age LMC clusters, at least 
for [Na/Fe], it should not be surprising if the range of the vari- 
ation is smaller than in the Galactic GCs (even setting aside 
the previous discussion of [O/Fe]). The trend of f p with mass 
found for Galactic GCs suggests that lower mass GCs are able 
to retain only a fraction of the AGB ejecta from first genera- 
tion stars. If these LMC clusters were not dramatically more 
massive in their past, then the expected range in the elemental 
abundance variations would be smaller, as the shallower po- 
tential wells would not be able to retain as much AGB ejecta. 
The point is that one should not expect clusters in other en- 
vironments to show the exact same abundance variations as 
observed in Galactic GCs. However, it is expected that any 
cluster showing an internal variation in age, as seen for these 
LMC clusters, should also show some variation in their light 
element abundances, as observed. 

Two important conclusions emerge from this section. First, 
statistically-significant star-to-star scatter in the [Na/Fe] abun- 
dances exists with in three of the four interm ediate-age LMC 
clusters studied bv lMucciarelli et all (120081) . This contradicts 
the conclusions of Mucciarelli et al., though they present no 
quantitative argument for a lack of star-to-star scatter. Sec- 
ond, a strong Na-O anti-correlation is not expected on theo- 
retical grounds for moderate-metallicity clusters because the 
expected range in [O/Fe] is small. For moderate-metallicity 



clusters, interpretation of any internal abundance variations 
must take this point into account, and it may therefore be more 
profitable to focus on obtaining [Na/Fe] abundances in such 
clusters. 

5. DISCUSSION 
5.1. Origins & Implications 

The preceding analysis has revealed several important facts 
related to the multiple stellar population phenomenon in GCs. 
Perhaps most important is the result that the fraction of 
present day GC mass comprised of pure AGB ejecta, f p , is 
large and strongly correlated with GC mass. The small scat- 
ter in f p at fixed mass suggests that whatever process shapes 
the elemental abundance trends is driven primarily by inter- 
nal processes, rather than external ones such as tides from the 
host galaxy. It is puzzling why f p is so strongly correlated 
with mass and yet such a weak function of mass (varying by 
less than a factor of two over two decades in GC mass). 

There are at least two possible explanations for the observed 
correlation between f p and GC mass. First, the amount of 
AGB material retained within a GC may decrease with de- 
creasing GC mass. This is plausible, as lower mass GCs have 
lower escape velocities, and yet the typical velocity of AGB 
ejecta is independent of GC mass. A second possibility is that 
GCs retain all of their AGB ejecta but this material is more 
strongly diluted by pristine material accreted from the ambi- 
ent ISM in lower mass GCs. Lower mass GCs could more ef- 
ficiently accrete pristine material if ISM s weeping is the dom- 
inant mechanism bringing in new material ( Conroy & Spergel 
120111) . The former explanation seems more natural and is 
therefore preferred. 

The derived fraction of mass in second generation stars 
presents one of the most significant outstanding puzzles. It 
is observed to be > 50% of the total present GC mass and 
is relatively independent of mass. This fraction is derived in 
the present work, but has been fo und also by many authors 
throug h different techniques (e.g., Smith 1987; Carret taet al.l 
1201 Pel) . If GCs were > 10 times larger at birth, then the early 
contribution of second generation stars to the total was small, 
of order a few percent. The fact that present day GCs al- 
ways end up with roughly the same fraction of second gen- 
eration stars strongly suggests that there is some mechanism 
that drives all clusters toward this state. One scenario that 
could achieve this is as follows. Imagine that some mecha- 
nism causes the first stellar generation to become unbound on 
the same timescale as the formation of the second generation 
(see below for an example). The second generation, being 
more tightly bound, would remain so, and might also be ca- 
pable of capturing of order its own mass in first generation 
stars. 

The second major outstanding puzzle is the more funda- 
mental issue of how and why the ancient GCs were initially 
so much more massive at birth. Standard dynamical effects 
cannot explain this, essentially because the relaxation time 
increases with increasing mass, and yet the mass enhance- 
ment factors are larger at higher masses. The only somewhat 
plausible scenario that has been proposed to explain this is 
primordial mass segregation of the s tars in the first stellar 
generation (e.g., D'Frco le et a I J 120081) . If a significant frac- 
tion of high mass stars were born near the cluster center, 
then as they evolve and die they will carry away a signifi- 
cant fraction of the total binding energy of the young GC, 
causing it to expand significantly. There is some circum- 
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stantial evidence for primordial mass segregation in nearby 
young clusters whose ages are comparable to their cross- 
ing times (e.g., [H illenbra nd & Hart mann 1998; d e Griis et al.l 
2003 lStolte et alJl2006tlGennaro et al.ll201 ll) . 

If primordial mass segregation were strong enough, it could 
unbin d young GCs on a timescale of < 1 Gyr (Vesperi ni et alJ 
2009), or at least le ad to a substantial amo unt mass loss and 
cluster expansion (Mar ks & Kroupal LTOlO). If this were a 
common feature of GC evolution, then perhaps GCs owe their 
very surviv al to the formation of tightly bound second gener- 
ation stars (D'Antona & Ventura 2008). These second gener- 
ation stars, born at the bottom of a deep potential well set by 
the first generation stars, should be relatively immune to the 
effects of primordial mass segregation. Detailed simulations 
will be req uired to validate this scenario, of which the simu- 
lations by (D'Ercole et al. 2008J) are an important first step. If 
confirmed, this would constitute a major revision to our un- 
derstanding of the survival of massive star clusters. 

A novel implication of the fact that the ancient GCs were 
much more massive at birth is their potential contribution 
to the Galactic stellar halo ( see also IVesperini et al.1 120101; 
ISchaerer & Charbonnell 1201 ll) . There are currently ~ 150 
known GCs (lHarri si ll 996). and they have an average present 
mass of 10 5 5 Mq. If these GCs had on average 10 times 
more long-lived stars at birth, then in total they would have 
equaled a mass of pa 5 x 10 8 Mq. This value is within a factor 
of a few of t he estimated stellar mass of the Galactic stel- 
lar halo (e .g.. ISiegel et al.l 120021) . and of the stellar halo in 
M31 (e.g., Kali rai et al.l l2006). This is the contribution to the 
stellar halo of GCs that remain bound at the present epoch. 
iMartell & Grebell (120101) have shown that approximately half 
of the Milky Way (MW) stellar halo could have formed from 
GCs that are now disrupted. An understanding of the early 
evolution of GCs may therefore be required in order to un- 
derstand the hierarchical formation of the stellar halos around 
galaxies. 

If the average GC today was 20 times more massive at 
birth, then it would have had a mass of pa 6 x 10 6 Mq, and 
it would have formed out of a giant molecular cloud (GMC) 
with a mass of at least ~ 10 7 Mq, assuming £sf = 0.5. Such 
large GMCs should form in the ma ssive, gas-rich proto galac- 
tic disks common at high redshift (Escal a~& Larsonl2008l) . In- 
deed, observations of disk-dominated galaxies at z ~ 2 have 
identified large numbers of super-star forming clumps with 
masses of order 10 9 Mq dGenzel et al.ll2.Q08h . Each of these 
clumps could easily spawn several massive young GCs. Even 
at the present epoch young clusters have been found with dy- 
namical masses in exces s of 10 7 Mq (e.g., iMaraston et alJ 
12004 iBastianet ail 12006). While rare at the present epoch, 
the conditions at high redshift may well have favored the for- 
mation of many such massive objects. 

M31 contains a large number of extended, low surface 
brigh tness GCs with half-light radii > 10 pc dHuxor et alJ 
2005) that have no counterparts in the Galaxy. The origin 
of these extended GCs is not known. In the context of the 
present discussion, they could arise from the rapid expansion 
caused by primordial mass segregation. If for some reason 
these clusters formed in a much more benign tidal field than 
the GCs in the Galaxy, then the loosely bound stars may not 
have been stripped by the present epoch. If this scenario is 
correct, then these extended GCs are comprised principally 
of first generation stars, and should contain a small core of 
second generation stars at their center. Radial gradients of the 



CN absorption feature in these GCs would be very interesting. 

The conclusion that GCs were much more massive 
in the past may also require some revision to dy- 
namical mo dels for the long-term e volution of the GC 
population (iGnedin & Ostrikerl 119971: iFall & Zhangl 120011: 
iMarks & Kroupall2010l) . 

While the present analysis provides a self-consistent ex- 
planation for the variation in the extent of the Na-O anti- 
correlation from cluster to cluster, complications arise when 
considering other light elements. In particular, some clus- 
ters that show an Na-0 anti-correla tion do not harbor a corre- 
sponding Mg-Al anti-correlation (Carrettaet al. 2009b). For 
example, the clusters NGC 6121 and NGC 6752 have very 
similar [Fe/H], total mass, and f p parameters, and yet the for- 
mer shows no star-to-star variation in either Mg or Al, while 
the latter shows clear variations in both. The Mg and Al yields 
depend on AGB mas s in a manner di fferent from the Na and O 
yields (IVentura & D' Antona 2008b), and so it is possible that 
differences in Mg-Al anti-correlations in GCs that contain the 
same Na-0 anti-correlation arise from a difference in the av- 
erage AGB mass t hat don ates ejecta for later star formation. 
See lCarretta et alJ (l2009bl) for further discussion. 

5.2. Alternative Explanations 

A number of alternative explanations have been proposed 
for various aspects of the standard model considered herein. 
Several of these will be briefly discu ssed in this section; the 
reader is referred to iRenzinil (120081) and IConroy & S pergel 
(1201 ll) for further discussion. 

While AGB ejecta is the favored source for the processed 
material, other proposed so urces include winds fro m massive 
(> 20 Mq) rotating stars (Decres sin et al.l l2007bl) and mas- 
sive binary star interactions dde Mink et alj |2009'). The most 
serious objection to these scenarios is that they are associated 
with massive stars with short lifetimes. It is therefore not nat- 
ural for the processed material to remain free of supernovae 
contamination. For example, it is not obvious why the second 
stellar generation should have the same [Ca/Fe] abundances 
as the first gene ration in these scena rios, and yet this is ob- 
served to be so (Ca rretta et al]l2010bT) . 

The similar 10 8 yr timescales for both intermediate-mass 
AGB stars and the drop in UV photons from the first gen- 
eration provide a simple framework to understand not only 
the ancient GCs but also the observed age spread within 
intermediate-age LMC clusters. Invoking processed material 
from massive stars requires either a more nuanced or an alto- 
gether different mechanism to be at work in the LMC clusters 
compared to ancient GCs. This seems unappealing on the 
principle of simplicity. 

Several authors have considered the possibility that stars 
with anomalous abundance ratios simply had their sur- 
faces polluted by AGB ejecta (e.g., iD'Antona et alJ 11983b 
iThoul et ani2002l) . In this scenario there is only one genera- 
tion of star formation and the observed range in light element 
abundances is due to the differing amount of surface pollution 
from star to star. A major advantage of this scenario is that the 
GCs need not have been substantially more massive at birth 
because a much smaller amount of AGB ejecta is needed to 
cover the surfaces of stars. 

This scenario is almost certainly ruled out by the lack of 
strong [O/Fe] and [Na/Fe] abundance variations between the 
main sequence and red giant branch (RGB). A 1 Mq star 
on the RGB has a convective envelope comprising approxi- 
mately 40% of its mass, while the same star on the main se- 
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quence has a negligible convective envelope. Therefore, this 
model would predict that as such a star evolves from the main 
sequence onto the RGB, surface pollution would be heav- 
ily diluted as the convective envelope deepens. This is not 
observed, which either m eans that the scena rio is to be dis- 
counted (as suggested in iCohen et al.1 12002). or the amount 
of surface pollution was substantial, so that it could not be 
diluted even by deep convection. The latter possibility then 
requires large amount of AGB ejecta, which means that large 
mass enhan cement factors are require d. 

Recently. iConrov & Spergell ( l201 11) proposed a scenario for 
the formation of multiple stellar populations within GCs that 
contains many of the ingredients listed in the Introduction, 
e.g., AGB ejecta and matter accreted from the ambient ISM 
as the source of second generati on stars. The pr i nciple differ- 
ence between that model and the D'Erco le et al.1 ([2008) model 
is that the former explicitly attempted to avoid the conclusion 
that the ancient GCs were much more massive at birth. In- 
stead, the large amount of material needed to form the sec- 
ond generation stars came primarily from the ambient ISM. 
As found in previou s work and confirmed herein (see also 
iD'Ercole et al.l20iTl) . the difficulty with such a scenario in ex- 
plaining the properties of the ancient MW GCs is that too few 
anomalous stars are produced. This is because the accreted 
ambient ISM, which has normal abundance patterns, domi- 
nates the mass budget. This does not preclude the possibility 
that star clusters in other systems formed a second generation 
of stars from primarily accreted ambient ISM material. 

Finally, several related scenarios appeal to highly unusual 
stellar configurations in o rder to expla i n the high number 
of anomalous stars (e.g., iBekkil 120061; iBekki et al.1 120071; 
iMarcolini et alj|2~009l) . These scenarios have two important 
properties in common. First, due to the special configu- 
rations required, they cannot be expected to operate at the 
present epoch, and therefore the multiple stellar population 
phenomenon observed in intermediate-age LMC clusters re- 
quires a second, distinct explanation. Second, they do not 
require present day GCs to have been substantially more mas- 
sive at birth. 

The scenario outlined by IBekkil (120061) . for example, en- 
visions GC formation at the centers of their own dark mat- 
ter halos. Stars throughout the dark halo are allowed to do- 
nate their AGB ejecta to the central regions where the GC 
is expected to form. Unfortunately, a number of indirect ar- 
guments disfavor G C formation at the center s of their own 
dark halos (see, e.g.. lConroy & Spergel l201 ll for a review), 
and direct dynamical searches for dark halos around present 
day isolated GCs strongly disfavor the prese nce of dark halos 
dBaumgardt et alj|2009t IConrov et alj|201 ll) . Since the deep 
potential well provided by a dark halo should efficiently re- 
tain supernovae ejecta, it is also not clear why GCs in this 
scenario should not all have large internal variation in heavier 
elements such as Fe and Ca. 

6. SUMMARY 

A simple model has been presented to interpret the ob- 
served Na-O anti-correlation in 20 Galactic GCs. The model 
assumes that GCs are composed of two generations of stars: 
a first generation formed from material with abundance ratios 
characteristic of field stars, and a second generation formed 
from a mix of AGB ejecta from the first generation and ma- 
terial accreted from the ambient ISM. Principal results from 
this analysis include the following: 



• The fraction of present GC stellar mass comprised of 
AGB ejecta is strongly correlated with GC mass, vary- 
ing from 0.2 to 0.45 over a factor of 100 in GC mass. 
This result is grounded in the observation that the ex- 
tent of the observed Na-O anti-correlation is strongly 
correlated with GC mass. 

• The fraction of GC mass in pure AGB ejecta, in con- 
junction with several well-motivated assumptions, pro- 
vides strong constraints on the composition and mass- 
loss history of Galactic GCs. The fraction of mass con- 
tained in second generation stars is always greater than 
50%, for GCs ranging in mass between 10 4 5 M and 
10 6 ' 5 Mq. The population of first generation stars in 
GCs must have been factors of 20-60 larger at birth 
compared to the present epoch. However, owing to 
the fact that first generation stars are subdominant at 
the present epoch, the total GC mass was only fac- 
tors of 10-20 larger at birth. These factors are mass- 
dependent and lower limits. The ancient GCs were 
therefore much more massive at birth. 

• Elemental abundance data on four intermediate-age 
(1-2 Gyr) LMC clusters are reinterpreted in the context 
of the models presented herein. It is found that three of 
the four clusters show unambiguous evidence for inter- 
nal variation in [Na/Fe] abundances. [O/Fe] values do 
not show signs of star-to-star variation, but this is shown 
to be a natural expectation in moderate metallicity clus- 
ters. The scenario invoked to explain the properties of 
the ancient Galactic GCs therefore appears to also be at 
work at the present epoch in the LMC. 
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